Survival of the d-wave superconducting state near the edge of antiferromagnetism in 

the cuprate phase diagram 
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In the cuprate superconductor YBa2Cu306+ x , hole doping in the Cu02 layers is controlled by both 
oxygen content and the degree of oxygen-ordering. At the composition YBa2Cu3 06.35, the ordering 
can occur at room temperature, thereby tuning the hole doping so that the superconducting critical 
temperature gradually rises from zero to 20 K. Here we exploit this to study the c-axis penetration 
depth as a function of temperature and doping. The temperature dependence shows the d-wave 
superconductor surviving to very low doping, with no sign of another ordered phase interfering with 
the nodal quasiparticles. The only apparent doping dependence is a smooth decline of superfluid 
density as T c decreases. 

PACS numbers: 74.72.Bk, 74.62. Dh, 74.25.Nf, 74.62.-c 



The high temperature superconductivity puzzle has 
important clues in the form of ordered states of mat- 
ter. One is the Mott insulator, where strong repulsion 
between electrons fixes them to lattice sites, with spins 
ordered antiferromagnetically (AF). When doped with 
holes, a superconductor is encountered, one in which 
paired holes condense into a superfluid, but bear the 
stamp of the parent Mott insulator by exhibiting a d- 
wave symmetry (dSC) favoured by strong electronic re- 
pulsion. This potential for finding new states of matter 
drives much of the interest in systems with strong elec- 
tronic correlations and in the cuprates much of this focus 
lies on the border between the Mott insulator and su- 
perconductor. In this regime of strong correlations and 
low hole-doping, other exotic phases have been predicted 
HI H IE [J IE j but experimentalists are hampered 
by a lack of homogeneous samples in which doping spans 
this range. Here we exploit a breakthrough [a| in the 
growth of high purity crystals of YBa2Cu306+ x , where 
doping can be continuously tuned with room temperature 
annealing. In the experiment presented here, the c-axis 
penetration depth A c is measured on a sample in which 
critical temperature T c is tuned from 4-20 K. Through- 
out this range, the temperature dependence A C (T) shows 
that the nodal quasiparticles of the d-wave superconduc- 
tor survive to very low doping, with no sign of another or- 
dered phase interfering with their behaviour. The change 
that does occur with doping is a smooth decline of su- 
perfluid density as the T c is tuned down to 4 K. 

A major advantage of working with the YBa2Cu306+ x 
system is that the hole-doping can be reversibly changed 
by controlling oxygen in a CuO chain layer situated 
0.42 nm away from the CUO2 planes that are the seat of 
the interesting physics. When the chains are completely 
depleted of oxygen, at YBa2Cu306, the Mott insulator 
is encountered and when they are filled YBa2Cu307 is a 
d-wave superconductor with a critical temperature (T c ) 
near 90 K. Between these extremes, Veal et al. || showed 
that the loosely-held dopant oxygens are mobile at room 



temperature and their gradual ordering into chain struc- 
tures pulls electrons from the CuC>2 planes ^J], increas- 
ing hole-doping and T c over time. An extreme version of 
this effect occurs in a narrow window of oxygen content 
near YBa2Cu30g.35, where crystals quenched from over 
570 °C initially do not superconduct, but become super- 
conducting with a T c that increases with time if they are 
allowed to order at room temperature. This opens the 
door to experiments over a wide range of hole-doping, all 
on the same crystal, with no change in cation disorder. 

Here we have exploited this in a measurement of the 
c-axis penetration depth, which provides information on 
the temperature and doping dependence of the superfluid 
density. A key feature of the c£-wave superconducting 
state is that the superfluid density has a linear temper- 
ature dependence due to the presence of nodes in the d- 
wave energy gap, in contrast to the exponential temper- 
ature dependence coming from the uniform gap of a con- 
ventional superconductor. This linear term is seen in mi- 
crowave measurements of the magnetic penetration depth 
X a b for superfluid screening currents flowing in the Cu02 
planes, which provide a direct measure of the superfluid 
density [Jjj . A common measurement geometry involves 
a microwave magnetic field applied parcllcl to the ab- 
plane of a crystal that is thin in the c-direction, thus 
avoiding large demagnetizing effects. However, this ge- 
ometry includes an admixture of A c , the magnetic pene- 
tration depth for current-flow perpendicular to the Cu02 
planes, which becomes dominant when A e becomes ex- 
tremely large at very low hole doping. So, here we turn 
our focus to measuring A c at low doping, since 1/A^ is 
related to superfluid density, although it is slightly less 
direct because it also depends on the coupling mechanism 
for currents between the Cu02 planes. 

Measurements of A c were made in a 22.7 GHz cylin- 
drical superconducting cavity |l2T | with the microwave 
magnetic field applied along the a-axis of a slab that 
was cut from a high purity crystal of YBa2Cu30g+ x 
[lflj . The sample was cut and polished to dimensions 
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(a x b x c) = (1.803 x 0.203 x 0.391mm 3 ) that make 
contributions from A Q f, negligible. The temperature de- 
pendence of the c-axis penetration depth AA C (T) = 
A C (T) — X C (T ) is extracted by measuring the shift in the 
cavity's resonant frequency as the sample temperature is 
increased above a base temperature T Q near 1.2 K. Im- 
portantly, A c is large enough that the microwave fields 
penetrate a substantial fraction of the crystal when it is 
in the superconducting state, but completely penetrate 
the sample above T c , making it possible to determine 
A C (T ) ~ A c (0) from the frequency shift seen when the 
sample is warmed through T c . A standard cavity pertur- 
bation relation (u> — lj)/u> = T[l — 26/dt&nh(d/25)], is 
used to extract the effective skin depth 6 from the com- 
plex resonant frequencies to and lo d of the cavity with 
and without the sample, respectively. T is a geometri- 
cal constant obtained by an independent measurement 
on a Pb-Sn sample cut to the same dimensions as the 
crystal and d is the thickness. The superfiuid screen- 
ing length is extracted from the effective skin depth via 
A~ 2 = Re{<5 -2 } + Lu 2 e r /c 2 , where the term containing 
the c-axis dielectric constant e r is a small contribution 
that displacement currents make to the screening. The 
value e r = 20 determined from far infrared measurements 
|l4j | comes mainly from phonons and at a frequency of 
22.7 GHz this gives a displacement current contibution 
of 4.5xl0 6 m~ 2 , a correction of order 1% except for the 
very lowest T c samples and very near T c . A weak tem- 
perature dependence of V due to demagnetizing effects 
has been corrected for. 

Fig. ^ shows the rise of both T c and 1/A 2 as the hole 
doping increases. The first remarkable feature of the data 
is the wide range of doping that can be accessed in just 
one crystal, with T c 's ranging from 4 to 20 K, all with 
sharp transitions of width 1 K. The three sets of curves 
indicated by blue squares, green triangles and red circles, 
show the results of three different stages in which the oxy- 
gen content was adjusted in a narrow range near O6.35 by 
annealing between 894 and 912 °C in flowing oxygen, fol- 
lowed by a homogenization anneal at 570 °C in a sealed 
quartz ampoule with ceramic at the same oxygen content, 
then a quench to °C. At each oxygen content the hole 
doping drifted up by way of gradual oxygen ordering at 
room temperature, while remaining mounted in the mi- 
crowave apparatus. A noteworthy feature of the data set 
is that the results at slightly different oxygen contents 
overlap one another over the course of the chain-ordering 
stage, indicating that for any particular T c , the behaviour 
of 1/A 2 (T) is dependent only on the hole-doping level, 
not on the details of the degree of chain order and oxy- 
gen content. A second striking feature is that at low 
temperatures the absolute change in 1/ A 2 , with tempera- 
ture has almost no doping dependence, in contrast to the 
large changes in 1/A 2 (0) and T c . Fig. Elfocuses on this by 
plotting 1/ A 2 (T) — 1/ A 2 (0) , which reflects the depletion of 
the superfiuid density when quasiparticles are thermally 
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FIG. 1: (colour online) Growth of critical temperature and c- 
axis screening with increasing hole doping in YBa2Cu306+ x , 
x~0.35. 1/Ac depends on the density of superconducting 
charge carriers that participate in the screening of applied 
fields and has been determined by a cavity perturbation mea- 
surement at microwave frequencies. The three sets of curves 
indicated by squares, triangles and circles, correspond to 3 
slightly different oxygen contents and at each oxygen content 
the hole doping was allowed to drift up by way of gradual 
oxygen ordering at room temperature. 



excited out of the superconducting condensate. There is 
a striking collapse of the data to a power-law temperature 
dependence at low temperature, with very little doping 
dependence. 

Qualitatively, the first conclusion that can be drawn 
from the power law behaviour is that the nodes charac- 
teristic of a d-wave energy gap continue to govern the 
low temperature properties. If a phase transition to an- 
other state with a gap opening at the nodes had occurred, 
one would see the appearance of exponential rather than 
power-law behaviour at low temperatures. The lack of 
such a change rules out transitions to other exotic super- 
conducting states such as d + is or d + id, which break 
time-reversal symmetry and open gaps at the nodes. The 
power law is close to T 2 5 , the same as that seen in mea- 
surements of A C (T) in YBa2Cu30g+x at much higher hole 
doping [J^- The absolute magnitude of the T 2 5 term is 
remarkably uniform, indicating that the nodal quasipar- 
ticles that control the low temperature properties of the 
d-wetve state are unaffected by any other physics devel- 
oping as the material approaches the limit of supercon- 
ductivity in the phase diagram. 

The T 2 - 5 power law for 1/A 2 (T) is different from the 
linear temperature dependence seen in A 2 fc (T). This 
near-quadratic power law has been seen in measure- 
ments of c-axis penetration depth in La 2 _ x Sr x Cu04+,s 
[13, Bi 2 Sr 2 CaCu 2 8 0, HgBa 2 Ca 2 Cu 3 8+5 H3, as 
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FIG. 2: (colour online) A plot of l/Aj?(T) - 1/A?(0) focusses 
attention on the depletion of the superfluid screening from its 
low temperature value, reflecting the loss of superfluid density 
as quasiparticles are thermally excited out of the supercon- 
ducting condensate. The data for all of the doping levels 
shown in Fig. nearly collapse onto a single curve at low T 
with a power law close to T 2 ' 5 . The inset shows the same 
data over a wider temperature range that spans the highest 
T c . The solid curve is a fit to the low temperature data of 
the higher T c samples, with a model that involves nodal quasi- 
particles in a d-wave superconductor, scattered as they tunnel 
from one CuCh plane to the next. 



well as in YBa2CusOfi-i-y at higher oxygen dopings of x 
= 0.60, 0.95 and 0.99 [15|. Understanding this difference 
from the afe-plane behaviour requires a description of how 
holes tunnel between CuC>2 planes. When constructing 
such a model, it is important to note that if the tun- 
neling preserves the momentum of the hole and is also 
momentum- independent, then one obtains for 1/A 2 (T) 
the linear T dependence that is seen in l/A 2 b (T). In- 
cluding scattering in the c-axis tunneling is a natural 
way to obtain a non-linear T-dependence, especially in 
light of the very high c-axis resistivity in these materials. 
The near-quadratic power law has been modeled with 
an anisotropic scattering mechanism for c-axis transport 
[Tg| . but recently Sheehy et al. have arrived at a detailed 
microscopic model for such a process [2(J. Within this 
model, the T-dependence depends crucially on the nodal 
quasiparticles that are a defining feature of the d-wave 
superconducting state. Written in terms of vector com- 
ponents relative to the nodes, these quasiparticles have 
an energy spectrum given by E p = y/ppv F + p A v A , a 
Dirac spectrum where pf and pa denote components of 
the quasiparticlc momentum perpendicular and parallel 
to the Fermi surface, vf is the Fermi velocity, and the 
gap scale va = (dA/d<j>)/pF is the slope of the d-wave 
energy gap A(0) = A cos (20) as it increases with angle 
<f> away from the node. This Dirac spectrum takes the 



form of a very flattened cone, since vf is known to be 
much greater than va HH]. 

Scattering is introduced as random spatial variations 
in the tunneling matrix clement. Wc take the average, 
squared matrix element to be a Gaussian 

<l^ P | 2 H^exp[-(fc-p) 2 /A 2 ] (1) 

with overall magnitude t\ and width A that takes on 
small values so that there is only a slight change in mo- 
mentum when quasiparticles tunnel from plane to plane. 
With this assumption, the T-dependent part of 1/A 2 is 
calculated within standard BCS theory for a d-wave su- 
perconductor. The solid curve in Fig. 2 shows that the 
model describes the temperature dependence at low tem- 
peratures extremely well, with a single set of parame- 
ters for the scattering (h/A = 12 nm, t± = 26 meV) 
and parameters for the Dirac cone (vp = 3 x 10 5 m/s, 
Vf/va = 6.8). A doping-independent vf has been sug- 
gested by angle-resolved photoemission (ARPES) mea- 
surements 22], so the doping-independent Vf/va ob- 
served here might imply that va does not change sub- 
stantially through the low doping range. The key puz- 
zle that has been solved is that the nearly quadratic 
temperature dependence is crossover behaviour in the 
range uaA < ksT < ufA. At lower temperatures 
a higher power law dominates and at higher temper- 
atures it is linear in T, but if vf/va >> 1 there 
is a wide range where the temperature dependence is 
roughly quadratic |2(J. This model does not explain the 
much higher T 5 power law observed in the tetragonal 
material HgBa2Cu04 +( 5 (which has been attributed to 
momentum-dependent hopping) E3 , but it does explain 
the power law seen forl/A? in the majority of orthorhom- 
bic cuprates [II El Gil 

This conventional treatment of the temperature depen- 
dence of 1/A 2 in terms nodal quasiparticles in a d-wave 
superconductor completely fails to describe the doping 
dependence of 1/A 2 , which varies nearly quadratically 
with the sample's T c , as shown in Fig. [3J Lee and Wen 
[2l| have noted that the superfluid density at higher 
doping displays a paradoxical disconnection between the 
slope of the linear temperature dependence , w hose dop- 
ing dependence is thought to be weak [lj, |23j , and the 
zero-temperature value, which declines rapidly as the 
hole doping is reduced. The paradox lies in the fact 
that the zero-temperature value seems to count the small 
number of holes added to the Mott insulator, but the 
temperature dependence is consistent with a large d-wave 
gap on the large Fermi surface expected for a high density 
of electrons. Here at much lower doping this phenomenol- 
ogy is reflected in the doping-independent temperature 
dependence of Fig. [21 and the dramatic doping depen- 
dence of 1/A 2 (0) in Fig. |3J Following a suggestion by 
Ioffe and Millis |26J, this disconnection can be modelled 
with an effective charge for the quasiparticles that falls 
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to almost nothing as the magnetism is approached. 
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FIG. 3: (colour online) Relationship between T c and the zero- 
temperature value of the c-axis screening length. The experi- 
mental correlation between these two quantities (circles) can 
be explained with the model of nodal quasiparticles cut-off 
by shrinking Fermi arcs. The inset shows the linear relation- 
ship between the cut-off energy E c needed to model l/A^ and 
the measured critical temperature T c . This model produces 
a nearly quadratic dependence between and T c , as indi- 
cated by the solid curve in the main figure. The few degree 
discrepancy between the T c predicted by the model and the 
measured T c is likely a consequence of critical fluctuations 
that govern the behaviour near T c . 



to zero for states away from the vicinity of the d-wave 
nodes, a view inspired by ARPES measurements that 
show the Fermi surface restricted to small arcs of well- 
defined quasiparticles when doping is decreased |24l |25| . 
The superfluid density then comes from a truncated por- 
tion of Fermi surface |2^| and Sheehy et al. [2(j have 
suggested that the doping dependence can be modeled 
with a Dirac cone that is cut off at progressively lower 
energy E c as doping is decreased. The inset in Fig. |2]dis- 
plays a linearly falling E c that is extracted from 1/Aj?(0) 
as doping and T c decrease, indicating a superfluid den- 
sity that steadily falls towards zero as T c does. What is 
remarkable is that in the face of this depletion of super- 
fluid density, the d-wave quasiparticle excitations nearest 
to the nodes survive and continue to govern the temper- 
ature dependence. It has recently been suggested that 
these nodal quasiparticles, seen here in samples with T c 
as low as 4 K, even persist to dopings below the supercon- 
ducting phase 22] and show up in thermal conductivity 
measurements in magnetic fields where T c is driven to 
zero [2^ | . Thus the entire superconducting portion of the 
phase diagram, and perhaps even more, is governed by 
nodal quasiparticles that survive even though the Fermi 
surface is reduced to small arcs and the superfluid shrinks 
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